Remote atom entanglement in a fiber-connected three-atom system 
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An Ising-type atom-atom interaction is obtained in a fiber-connected three-atom system. The interaction is 
effective when A as y » g. The preparations of remote two-atom and three-atom entanglement governed by 
this interaction are discussed in specific parameters region. The overall two-atom entanglement is very small 
because of the existence of the third atom. However, the three-atom entanglement can reach a maximum very 
close to 1. 
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I. INTRODUCTION 



Generating the entanglement between spatially separated 
atoms plays an important role in quantum information 
processing and quantum computation, such as quantum 
storage' 11 , quantum key distribution' 21 and quantum states 
swapping 131 . To efficiently entangle two or more distant 
atoms, one must create some kind of direct or indirect interac- 
tion between them, such as by adopting appropriate measure- 
ment on optical fields that conditionally interact with atoms 
and thereby the atoms (as a subsystem) can be projected to 
an entangled state, or by using quantum-correlated fields in- 
teracting with atoms and thereby the entanglement among the 
fields can be transferred to atoms. Based upon this, a vari- 
ety of schemes for entangling distant atoms or distant photons 
have been proposed recently' 4-141 . For example, fascinating 
schemes have been presented to efficiently entangle distant 
atoms, where the single-photon interference effect was ap- 
plied with' 41 or without' 51 weak driving laser pulse. Recently, 
S. Mancini and S. Bose proposed a novel scheme to directly 
entangle two atoms trapped in distant cavities' 61 which were 
connected via optical fibers. Using input-output theory, un- 
der adiabatic approximation, the authors obtained an effective 
Ising model for two atoms. In their scheme, photon acted as 
an intermediate quantum information carrier and mapped the 
quantum information from the atom in one cavity to that in 
another. Such systems are meaningful not only in quantum 
measurement or testing Bell's inequalities but also in poten- 
tial applications such as quantum encryption' 151 or construct- 
ing universal quantum gates 1161 that are essential for design- 
ing quantum network. Nevertheless, in discussing quantum 
networking with trapped atoms and photons in cavity QED 
system' 171 , two problems should be overcome: How to gen- 
erate the entanglement of a N-atom system? What is the ex- 
act influence of the collective interaction on the entanglement 
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shared by remote atoms? These problems have been discussed 
intensively, for instance in the scheme proposed by Cabrillo 
et al' 4 '. The simplest multi-atom case is a three-atom sys- 
tem which might be an intuitive extension from a two-atom 
case. In our scheme, We extend the model of two-atom cir- 
cumstance in Ref. [6] to three-atom which turns out to be a 
three-atom Ising model. Such an approach might be mean- 
ingful in discussing the above problems for multiple distant 
atoms. We firstly investigate the dependence of the effective 
Ising coupling coefficients on the atom-cavity detuning and 
cavity leakage. Then, we discuss the influence of an atom on 
the other two atoms entanglement properties. Furthermore, 
we study the characters of remote three-atom entanglement 
and the tangle between one atom and the rest two atoms. 



II. OPTICAL FIBERS CONNECTED THREE-ATOM 
SYSTEM 

The schematic setup for our system is shown in Fig. 1. 
Three identical two-level atoms 1, 2 and 3 are trapped in spa- 
tially distant cavities C\, C2 and C3 respectively. All the cav- 
ities are assumed to be single-sided ones. Three off-resonant 
external driving field e\, £2 and £3 are applied upon cavity C\, 
C2 and C3 respectively. In each cavity, a local laser field that 
is resonantly coupled to the atom is applied. Two neighboring 
cavities are connected via optical fibers. Apparently, the sub- 
system constituted by cavities C\ and C2 or C2 and C3 is just 
the setup proposed in Ref. [6]. 

In the interaction picture, using cavity input-output 
theory' 18 ' and taking adiabatic approximation' 19 ', we obtain 
an effective Hamiltonian for this system as (see Appendix A) 

H e /f = Jnv\o\ + Jzso\o% + + ^ F,(o-+ + <r~), (1) 

i 

which is a three-particle Ising chain with magnetic fields 
perpendicular to the z direction' 201 , J\2 and ^23 represent 
the nearest-neighbor (NN) atoms coupling coefficients, while 
731 represents next-nearest-neighbor (NNN) atoms interaction 
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FIG. 1: Schematic setup of the system. Three two-level atoms locate in 
spatially separated single-sided optical cavities that are connected via optical 
fibers. 



strength. cr\(i = 1, 2, 3) is spin operator of atom i, crf(crv) is 
atomic raising (lowering) operators. F, is the magnitude of the 
locally applied laser field interacting with atom i. We define 



hi =2yoX 2 Im {a x a\^ /[M 2 - W 2 ]} , 
J23 =2y Q x 2 Im [a 3 a* 2 e ilhl /[M 2 - W 2 ]} , 
hi =2yoX 2 Im [y a 3 a* 1 e i(lh>+ * a) /[M(M 2 - W 2 )]} , (2) 

where yo is the cavity leakage rate, \ - g 2 /A, 8 is me cou- 
pling strength between the atom and the cavity field in cav- 
ity Cj, A is the detuning between the atomic internal transi- 
tion and cavity field frequency, where, large detuning approx- 
imation has been assumed, i.e. A » g, and M = i A + yo, 
W 2 = y 2| e ;(fa+<An) +e /(te+fe)j. The phase factors <f> u (ij = 
1 2, 2 1 , 23 , 32) are caused from the photons transmission along 
optical fibers from cavity Cj to cavity C, [21] . Physically, they 
depend on the frequency of the photons and the distance be- 
tween cavities. And 



a 2 = 



CK3 = 



s\M l + s 2 My e">" 2 + y^e™" - sie !<P| ] 
M(M 2 - W 2 ) ' 
e 2 M + y (eie'' fcl + e 3 e ilf ' 23 ) 
M 2 - W 2 ' 
s 3 M 2 + e 2 My e i ^ 1 + yg[eie i%1 - e 3 e i% ] 



M(M 2 - W 2 ) 



(3) 



where ©i 3 = <p l2 + <p 2 3,®i = <t>23 +te,®3i = + 021,^3 = 
021 + 012- The global system is now determined by a series of 
independent parameters as s i , s 2 , £3 , A, yo , <p \ 2 , ^32 , F 1 , and 
T3. In next section, we discuss the optimal region of the pa- 
rameters for the preparation of remote atom entanglement. 




FIG. 2: The amount of NN atoms coupling coefficients(/i2(/23)) under dif- 
ferent detuning (A) and cavity leakage rate (yo)- All the variables are pre- 
sented in unit of ?. 



III. PARAMETERS SPACE DESCRIPTION OF ISING 
COUPLING COEFFICIENTS 

From Eqs. (3), the condition M 2 * W 2 leads to large Ising 
coupling coefficients. Fhis condition also keeps the validity of 
the adiabatic approximation in case of weak local laser fields, 
i.e. J\ 2 {J 2 3, J31) » F,, i = 1, 2, 3. We can further simplify the 
condition as 



hi + <Pn 



*23 + <P32 



To- 



(4) 



For simplicity, we assume s\ — s 2 = £3 = eo, Fi = T 2 — 
F3 = Fo. The parameters space can now be expressed in unit 
of g as (A/g,yo/g)g. In Fig. 2-3, we give the description 
of Ising coupling coefficients for NN and NNN atoms in the 
parameters space. Where we assume eo = 2g. We can see that 
the coupling coefficients for NN atoms as well as NNN atoms 
can be divided into two regions: (a) the region where A > yo, 
(b) the region where A < yo. In most area of the two regions, 
the coupling coefficients are very small, only in the regions 
just besides the line A = yo are they large enough so that the 
validity of the adiabatic approximation can be kept. 

In the following discussions, we will study two-atom entan- 
glement nature and three-atom entanglement properties based 
on the parameter space. 



IV. NEAREST-NEIGHBOR AND 
NEXT-NEAREST-NEIGHBOR REMOTE TWO-ATOM 
ENTANGLEMENT 

In this section, we discuss the nature of remote two-atom 
subsystem entanglement which is generated in our system. 
Wootters proposed a general measurement for the amount of 
two-qubit (noted as 1 and 2) entanglement. It is named as 
Concurrence [22] : 



C12 = max{0, A\ - A 2 - A3 - A4}, 



(5) 
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FIG. 3: The amount of NNN atoms coupling coefficients^ i ) under different 
detuning (A) and cavity leakage rate (yo). All the variables are presented in 
unit of g. 



0.10 




gt 



FIG. 5 : Descriptions are same as those in Fig. 4 but for To = 0.2g. Note that 
the whole time scale is half of that in Fig. 4. 
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FIG. 4: Entanglement of atom 1 and 2 versus time for A: (solid line) /31 
l.2g, B: (dotted line) /31 x -l.2g. Where we assume To = O.lg. 



where /I, are the non-negative square roots of the four eigen- 
values of non-Hermitian matrix pnpn with defined as 
(cr y ® cr v ,)p* 2 (cr > , ® (T y ), where pu is the density matrix of the 
two-qubit system. 

We depict the two-atom entanglement situation in Fig. 4-5 
for different parameter spaces. We assume that all the atoms 
are initially in their ground state, so that |(A(0)) = \g)i ® \g)2 ® 

l*>3- 

Firstly, we investigate the influence of 731 on the entan- 
glement of NN atoms. In Fig. 4, we adopt appropriate val- 
ues of A and yo (which satisfy the condition in Eq. (4)) 
and assume To = Q.lg. The Ising coupling coefficients are 
J12 = J23 ~ -2Ag, /31 ~ l.2g for solid line, and /31 « -1.2g 
for dotted line. If all the signs of the coefficients are reversed, 
the resulting concurrences are not changed. Evidently, rela- 
tive larger entanglement for NN atoms can be obtained when 
JniJii) ■ J31 < 0. While, compared with the result in Ref. [6], 
the overall entanglement is very weak since two-atom subsys- 
tem is in mixed state during the evolution. 
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FIG. 6: Entanglement of atom 1 and 3 versus time when the laser field ap- 
plied in cavity C2 is turn off, that is Ti = T3 = 0.3g, T2 = 0. 



In addition, the NN atoms entanglement can be manipu- 
lated through the alternating of the locally applied laser fields. 
In Fig. 5, we adopt the same parameters as those in Fig. 4 
but for Tq = Q.2g. Fig. 5 indicates that, the increase of To 
remarkably improves the NN atoms entanglement. The pe- 
riod is depressed, but the amount of entanglement is much en- 
hanced. The amount of entanglement for NNN atoms, under 
the parameters we assumed, is generally much weaker than 
NN atoms. To improve the entanglement for NNN atoms, 
The Ising coupling coefficient between NNN atoms must be 
enhanced. In Fig. 6, we depict the entanglement for NNN 
atoms. Correspondingly, Jn = J23 = -2Ag, J31 = l.2g. To 
modulate the entanglement, we let F2 = 0. Under this circum- 
stance, the entanglement for NNN atoms can compare with 
that for NN atoms (see Fig. 4). 
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FIG. 7: The three-atom entanglement Ci23(solid line), the tangle 
Ci(23)(dotted line), and the Concurrence Ci2(dashed line) versus gt, where 
J\l(hi) = -2Ag, J31 = -1.2g and To = 0.2g. To distinguish Ci(23) from 
C123, the line for Ci(23) is raised to 0.1 + Ci(23). 



V. THE REMOTE THREE-ATOM ENTANGLEMENT 
PROPERTIES 

The intrinsic three-partite entanglement which is widely 
used for measuring three-partite entanglement of pure states 
is defined as [23] 

Cl23 = Ci(23) - Cj 2 - Cj 3 , (6) 



VI. CONCLUSION 

We have obtained a three-atom Ising chain in cavity QED 
system by connecting three distant cavities via optical fibers. 
The Ising coupling coefficients are found to be large in the 
region where A « jo » g, which keeps the validity of the 
adiabatic approximation. We have discussed the generation of 
remote atom entanglement. The overall two-atom entangle- 
ment is very small because of the existence of the third atom. 
While, the NN atoms entanglement can be improved when the 
coupling coefficient of NNN atoms has a contrary sign with 
respect to that of NN atoms. The locally applied laser fields 
play an important role in modulating the entanglement quan- 
titatively and qualitatively not only for NN atoms but also for 
NNN atoms. Furthermore, we have studied the remote three- 
atom entanglement and the tangle. It is shown that three-atom 
entanglement, which has a much longer period than two-atom 
entanglement, can reach a maximum very close to 1 . 

In addition, it should be noted that the dissipation of the 
photon information along the fibers should be investigated, 
while, the dissipation can be included in the Ising coupling 
coefficients and act as a decaying exponential factor e~ vL , 
where v is the dissipation rate per meter, L is the total length 
of the fiber' 251 . The phase factors e'^' 2 and e' te are then re- 
placed by e i4lll ~ vLa an( j e 'fe-vL 23 j n f actj me dissipative ef- 
fect along fibers can be compensated by lowering the detun- 
ing A. One can obtain large Ising coupling coefficients by 
adopting the parameters in the regions just besides the line 
A as V2e- v < Ll2+L2 3) - ly - 



where Ci(23), which represents the tangle between a subsystem 
1 and the rest of the global system (denoted as (23)), is written 
as 

C 1(2 3) = ADetp, = 2(1 - Trp\). (7) 

In Fig. 7, we plot the remote three-atom entanglement 
Ci23(the solid line), the tangle Ci(23)(fhe dotted line), and the 
Concurrence Ci2(the dashed line) for 1^(0)) = \g)\®\g)i®\g)i- 
Where the corresponding parameters are same as those of the 
dotted line in Fig. 5. In fact, there is only very little difference 
between C123 and Ci(23). To distinguish Ci(23) from C123, the 
line of Ci(23) is raised to 0.1 + Ci(23). 

It has been pointed in last section that the Concurrence 
that represents the bipartite entanglement between atom 1 and 
atom 2 is very small. While, the tangle, which expresses the 
entanglement of atom 1 and the rest of the global system, and 
the remote three-atom entanglement can reach maximum val- 
ues almost 1 for intermediate values of gt. Under the condi- 
tion of strong Ising coupling coefficients (Jn(J23, ^31) » To), 
if we express the tangle between atom 1 and the rest of the 
global system as a sum of the remote three-atom entanglement 
and the bipartite entanglement, the three-atom entanglement 
will act as the largest contribution. It has been concluded that 
in a system of N spin-half particles, under the condition of 
strong Ising coupling coefficients, the iV-partite entanglement 
will be dominant' 241 . 
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Appendix A 

In the interaction picture, the Hamiltonian of the global sys- 
tem can be written as 

Hint = Hi n ,\ + Hi nt 2 + Hi n ,3 + i/;„ ? 4, (Al) 

where H int \ represents the effective interaction of atoms and 
cavity fields, H in a is the coupling between external driving 
fields and cavity fields , i/,„,3 represents the interaction of lo- 
cally applied laser fields and atoms, H in ,4 is the interaction 
of cavity fields and their environment which is described as 
a superposition of series of harmonic oscillators. Under the 
condition of large detuning, we have [26] 

Hintl =xYj A ! A i° Z i> < A2 > 
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where i= 1,2,3, A,(At) represent cavity fields annihilation (cre- 
ation) operators in cavities C, . And [181 



H int2 = Y j s i (At+A i ), 



(A3) 



To solve these equations explicitly, we firstly obtain the ex- 
pectation values of cavity field operators through 



d{A x ) d{A 2 ) d{A 3 ) 



dt 



dt 



dt 



= 0. 



(A10) 



H int3 =Y j r i (<ri+a-). 



(A4) 



We assume r, are weak enough so that the quantum adia- 
batic theory' 191 can be applied in the following calculations. 

+oo 

HintA = i j ' dco^ i K Ci [bc i (.cS)A+ + h.c], (A5) 



where bc;(oj), i = 1, 2, 3, are the annihilation operators of the 
harmonic oscillators with frequency co. kq. are the interaction 
strengths between cavity C, and the harmonic oscillators. The 
kinetic equations for cavity field operators turn out to be [18] 

Ai = -(/A + ixcr\ + -^-)Ai + yfycTAun + e i> 

Tc - 
A 2 = -(/A + ix(T 2 + — -)A 2 + ^JcMm + £2, 



- 2 

rc 

2 



A 3 = -(/A + i X (r\ + -^-)A 3 + yfn~A 3M + £ 3 , (A6) 



where y Cl = 2n[k Ci {.tS}] 2 (i = 1,2,3). If cavities C\ and C 2 
are connected via optical fibers (as shown in Fig. 1), so are 
cavities C 2 and C 3 , the input-output conditions should be in- 
cluded, so that [21] 



A, = -2|la, + VmAa^e**, 

A 2 = -^A 2 + y/yFAloute 1 * 2 ' 

+ Vy^A 3 , „,e ifa , 
M = - J -^A 3 + ^A Xout e^\ 



(A7) 



For simplicity, assuming the decay rates y Cl = yc 2 = 7d = 
jo and taking into account the usual boundary conditions [18] 



A UoM + A Un = VroAi, 



(A8) 



where i = 1 , 2, 3, we can rewrite the kinetic equations for cav- 
ity field operators as 

Ai = -MA\ - ixA\(r\ + VroA Un 

+ ^ 12 ( ro A 2 - VroA2,i„) + £i, 
A 2 = -MA 2 - ixA 2 o\ + y/yoA 2Jn 

+ J] e^iyoAj - yfftAj*) + e 2 , 

7=1,3 

A 3 = -MA 3 - tfcA 3 £7^ + VToA 3 ,;„ 

+ (y A 2 - VroA 2i( „) + e 3 • ( A9) 



The steady states for cavity fields in C\ , C 2 and C3 can be 
obtained as 



a 2 = 



a 3 = 



SiM 1 + e 2 My e ilp12 + yl[s 3 e i& " - sie i<b '] 

M(M 2 - W 2 ) 
s 2 M + ro(eie i021 +g 3 e ite ) 
M 2 - W 2 ' 
s 3 M 2 + E 2 My Q e^ 2 + y^e 10 ^ - e 3 e^] 

M(M 2 - W 2 ) 



, (All) 



where 13 = 4>i 2 +(p 23 , Oi = (f> 23 +(/> 32 ,®3i = 4> 32 +(j) 2 i,<& 3 = 

4>2l + <t>\2- 

Then, in the regime of strong cavity leakage and large de- 
tuning (which lead to yo, A >> x), the kinetic Eqs. (A9) are 
reformed as the following homogeneous linear equations: 



d\ = -Mai ~ ix a \°~\ + Vro«M 

+ e"* 12 (yofl2 - VT0 a 2,i«), 
d 2 = -Ma 2 - ixa 2 a-\ + yfyoa 2J , 

7=1,3 

a 3 = -Ma 3 - ixa 3 cr\ + y/y^a 3J , 
+ e ,fe (y «2 - Vro«2,m), 



(A12) 



where we have replaced field operators A, with a, + a, 
(i=l,2,3). In solving Eqs. (A12), one can adiabatically elim- 
inate the effect of vacuum input noise. The resulting cavity 
field operators are now represented by linear combinations of 
atomic spin operators 07 (i= 1,2,3). Substituting the resulting 
field operators into Eq. (AI), we get the effective Hamiltonian 
of the global system in the interaction picture as 



H e ff = J\io\o\ + Jno\o% + J3i(r z 3 <r\ 
+ £r,K +£ rr), 



(A13) 



where J\ 2 , J 23 and /31 are expressed by Eqs. (2). 

In deriving Eq. (A13), we neglect self-energy terms includ- 
ing cr^ and self-interaction terms including and (crp 3 
that do not change the initial system state. Also, we elimi- 
nate higher order terms that include x^&\v\v% since the cor- 
responding coupling coefficients are much weaker than J\ 2 , 
J 23 and 731. The typical difference between this Hamiltonian 
and that in Ref. [6] lies in the third term in Eq. (A13). 
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